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ABSTRACT 

We examine in detail how BBN theory is constrained, and what predic- 
tions it can make, when using only the most model-independent observational 
constraints. We avoid the uncertainties and model-dependencies that neces- 
sarily arise when solar neighborhood D and 3 He abundances are used to infer 
primordial D and 3 He via chemical and stellar evolution models. Instead, 
we use 4 He and 7 Li, thoroughly examining the effects of possible systematic 
errors in each. Via a likelihood analysis, we find near perfect agreement be- 
tween BBN theory and the most model-independent data. Given this agree- 
ment, we then assume the correctness of BBN to set limits on the single 
parameter of standard BBN, the baryon-to-photon ratio, and to predict the 
primordial D and 3 He abundances. We also repeat our analysis including re- 
cent measurements of D/H from quasar absorption systems and find that the 
near perfect agreement between theory and observation of the three isotopes, 
D, 4 He and 7 Li is maintained. These results have strong implications for the 



chemical and stellar evolution of the light elements, in particular for 3 He. 
In addition, our results (especially if the D/H measurements are confirmed) 
have implications for the stellar depletion of 7 Li. Finally, we set limits on 
the number N u of neutrino flavors, using an analysis which carefully and sys- 
tematically includes all available experimental constraints. The value N u = 
3.0 fits best with BBN and a 95% CL upper limit of N u < 4 is established. 
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1 Introduction 



Because of the central role big bang nucleosynthesis (BBN) plays in the stan- 
dard cosmology, it is crucial to understand how robust the BBN results are. 
Consequently there has recently been intense scrutiny of different possible 
sources of uncertainty in the BBN analysis. There has been attention to 
refining the theoretical calculations, particularly the 4 He yields [|T], |2j , as well 
as quantifying their uncertainties due to the errors in their input parameters 
[0, [I], |, |^, 0], and there has been examination of the observational proce- 
dures and results for possible systematic errors || |J. Finally, there has been 
a focus on the model-dependence that arises in making the link between the 
observables and the theory predictions; in this paper we explore in detail an 
analysis[|l(| that minimizes this model-dependence. 

It is useful to situate our present concerns in the context of the subject 
of BBN analysis, which consists of three conceptually independent features. 
(1) The theory itself rests upon first-principles calculations and/or detailed 
laboratory results; the only really free parameter to the theory calculation 
is the baryon-to-photon ratio r\ = riv/n^, and the output for each value of r\ 
consists of the primordial abundances of the light elements D, 3 He, 4 He, and 
7 Li. (We will consider throughout only homogeneous BBN, and furthermore 
we will first study the "standard" model with N v = 3; later we will relax the 
latter assumption and look at possible constraints on N u .) (2) The observa- 
tions pertinent to BBN are determinations of the light element abundances 
in various astrophysical settings; for the most part, these environments are 
at low or zero redshift, i.e., contemporary to our own, with the exception of 
QSO absorption line system (QSOALS) D/H measurements which we dis- 
cuss in more detail below. As the observable abundances are contemporary, 
one must invoke (3) the final facet of BBN analysis, namely some account 
of galactic chemical evolution that models the evolution of elemental and 
isotopic abundances from their primordial state to the present. In this pa- 
per we emphasize that the chemical evolution results needed for each of the 
light isotopes have very different degrees of model-dependence; consequently 
we will retain only the results which are the most model- independent, and 
pursue their implications, both for BBN and for chemical evolution. 

It is necessary to distinguish between the general implications chemical 
evolution gives and its detailed results. Generally, metal abundances (i.e., 
abundances of elements with Z > 6) increase with time, while the specific 
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histories require the detailed, quantitative modeling of abundance evolution. 
Turning to the light elements, one finds that to infer the primordial D and 
3 He from local and recent abundances (i.e., solar system and the ISM) re- 
quires detailed modeling. Of the two, D has the fewest complications and 
depends only on the chemical evolution model features and not at all on 
stellar evolutionary data, since D/H is entirely destroyed in the pre-main se- 
quence phases of stars. 3 He is subject to considerable uncertainties regarding 
not only chemical evolutionary parameters but also its stellar processing; the 
possible production in low mass stars and the possible destruction through 
mixing effects in the giant branch stars. On the other hand, to infer primor- 
dial abundances for 4 He and 7 Li requires no detailed modeling, not because 
their evolution is simpler, but because we may determine their evolution 
empirically. Specifically, we observe 4 He and 7 Li in systems with different, 
and in particular low metallicities. Since metals increase with time, we may 
trace the evolution and infer the primordial abundances by extrapolating the 
observed trends to zero metallicity. This relatively simple procedure is, how- 
ever, complicated by the very likely presence of systematic errors in the data, 
and so we will address in detail the possible instances of and distributions of 
systematic errors. 

Given that the inferred primordial D and 3 He abundances suffer from 
large chemical evolution uncertainties, there are several approaches one may 
take when trying to constrain BBN. The traditional strategy has been to use 
the chemical evolution results for primordial D and 3 He, and either trying to 



quantify the concomitant uncertainty this introduces [11, O, 13], or to use 



a generalized form of chemical evolution which is designed to be "generic" 
[|14] , [l], [15| . This approach has the particular appeal that, as it uses all of 
the light elements, it strongly overconstrains the one-parameter BBN theory. 
On the other hand, there remains a concern about the chemical evolution 
framework itself, which plays a central role. As no first-principles calculations 
exist for this, it is hard to be certain that the idealizations in the available 
models can capture all relevant aspects of the (unavailable) full solution. 
This procedure breaks down when the data begin to be refined enough to 
show inconsistencies between theory and observation and the conclusions for 
BBN become overly sensitive to highly uncertain modeling. 

An alternative means of addressing the chemical evolution uncertainties 
is to minimize the reliance on chemical evolution by using only its most 
general features fTOf. That is, one uses only the inferred primordial 4 He 
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and 7 Li abundances, rather than all four light element isotopes. Of course, 
one always has to use chemical evolution at some level, but for 4 He and 
7 Li we need only assume that metal abundances (specifically C, N, O, and 
Fe) increase with time. With only this assumption, and a prescription for 
including possible systematic errors, one may constrain BBN. Note that in 
the case of 7 Li, which is observed in old low metallicity stellar populations, 
we also include an uncertainty on the degree to which these stars did not 
deplete their 7 Li abundance. 

In this paper, we will first review the observational status of all four light 
element isotopes. In the cases of D and 3 He, we will also review briefly the 
arguments which limit their implementation in an analysis of BBN. Thus 
we asses the validity of BBN using the most model-independent data. We 
quantify the success of BBN using likelihood analysis, including theoretical 
uncertainties as well as the observational statistical and systematic errors. 
This analysis is discussed in detail, and attention is given to the significance 
of the amplitudes of the likelihood functions as indicators of the goodness of 
fit. With this analysis we find, using only 4 He and 7 Li, that BBN theory is 
consistent with the observations, and we determine the best values for 77. We 
also attempt to quantify the reliability of these predictions. 

Though we are privileging 4 He and Li, this is not to say that these 
abundances are known with absolute certainty. We take as a starting point 
the abundances directly from the observations. These observations are in 
remarkable agreement for a baryon-to-photon ratio which is rather low, 77 ~ 
1.8 x 10~ 10 . A similar range has been also suggested in [16]. We will also 



show the effect of modifying the primordial abundances. For example, it is 
sometimes argued that the lithium abundance observed in halo stars is not 
the primordial value but represents a depleted primordial abundance [[[]], 
We will show how our results change when the effects of 7 Li depletion are 
included. 

As we have noted, the traditional use of solar and ISM data on D/H 
requires the application of chemical evolution before one may use it to con- 
strain BBN; however, there is a real chance that the recent measurements 



of D/H in quasar absorption systems Jl9], |2C§ are in fact true measures of 



the primordial D/H abundance. Therefore, we will also show results of the 
likelihood analysis which includes 4 He, 7 Li and D. In fact we find that the 
agreement is remarkably good, suggesting perhaps that the putative D in 
quasar absorption line systems is real and very close to primordial levels. 
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Having determined the success of BBN, we may then use it as has long 
been done, as a powerful tool in particle astrophysics. On the astrophysics 
side, we predict the abundances of primordial D and 3 He on the basis of the 
4 He and 7 Li data. These results are (in addition to being available for testing 
against quasar absorption line systems) then useful as inputs in chemical 
evolution models. We also determine the potential for BBN to set limits 
to 7 Li depletion; these constraints become considerably stronger, when D is 
included in the BBN analysis. On the particle side, we let N v vary, and 
quantify the allowed extra relativistic degrees of freedom. As we will show, 
if we were to use BBN, to "predict" the value of N u we would find, again 
remarkably, that the best value is N v = 3.0, independent of whether or not 
we include D in the analyses. To wit, the uncertainty in this "prediction" 
depends primarily on the uncertainty in the 4 He observations, which at the 
the 2 u level translates to a BBN upper limit of about 4 on N v . 

The paper is organized as follows: The observational data on each of 
the four light element isotopes is discussed in section 2, as is the role of 
chemical evolution in deriving primordial abundances of each isotope. We will 
carefully describe our choices of the input abundances and their statistical 
and systematic uncertainties. In section 3, we describe our likelihood analysis 
method and present results for BBN consistency in section 4. Included in 
section 4, is our discussion of the potential for constraining 7 Li depletion 
in stars, and our analyses including D with an assumed primordial value 
taken from the QSOALS measurements p5| . In section 5, we discuss the 
implication of these results on particle astrophysics constraints, particularly, 
the value of N v . Finally, in section 6, we draw our final conclusions. 



2 Data 

The literature on light element abundance measurements is large, and is 



well-reviewed in, e.g., references M, 21, 22]. Here we will give only a brief 



summary of the current best abundances for each light isotope, with em- 
phasis on new developments, particularly the possible detection of D/H in 
quasar absorption line systems (QSOALS). We will also discuss the issue of 
systematic errors in the data, and the possible forms such errors might take. 
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2.1 Deuterium 



As discussed in §1, until recently D abundances were only available for the so- 
lar system and the local ISM. The solar value for D/H = (2.6 ±0.6±1.4) x 1(T 5 
(see [E3] for a recent discussion) has a considerable uncertainty. Individual 



ISM observations are much more accurate [24| but they appear to show dis- 
persion along different lines of sight through the ISM |2"5| , PB"| . This result 
is unexpected in standard models of chemical evolution, and again suggests 
their incompleteness. Moreover, neither of these measurements is directly 
attributable to a primordial value. The D (and 3 He) in these systems has 
suffered a considerable degree of processing, and so determining a primordial 
abundance from these requires detailed chemical evolution modeling. As we 
have argued, our approach is to avoid such modeling, and so we will not 
employ these results. 

The observational situation for D is rapidly changing however, and is 
sparking a renewed interest in D evolution. Observers have employed high- 
resolution spectra of quasar absorption line systems to determine the D/H 
abundance via its isotope-shifted Lyman-a line. Initial reports |H| gave a 
high abundance, D/H ~ 2 x 10~ 4 , but there remained a significant chance 
of D being mimicked by an interloping H cloud, which would lead to an 



overestimate of the D abundance. Recently, new observations |2(J of this 



same system have in fact resolved this original system into two components 
each with an abundance comparable to the old result. The weighted average 
of these two measurements is [pOR , 



=(1.9±0.4)xl0- 4 (1) 

Vri/ 2=3.32 

The resolution into two components dramatically reduces the likelihood of 
an interloper, strongly suggesting that the observations are really measuring 
D. Caution is still warranted however, (as it is with any new observation) 
until the systematic uncertainty in these observations is clarified and more 
importantly until these results can be confirmed in the QSOALS in other 
directions. For these reasons, we do not include the D/H results in our first 
BBN analysis. We do however, repeat our analyses using the primordial 
value for D/H from Eq. ([!]). 
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2.2 Helium-3 



Abundances of 3 He are only available for the solar system and the ISM. The 
solar value for 3 He/H is (1.5 ±0.2 ± 0.3) x 10~ 5 and varies between (1 to 
5) xlO -5 in galactic HII regions |37| (see also [23] for a general discussion). 
Thus as in the case of D, one must use detailed chemical evolution models to 
extrapolate a primordial abundance. This procedure is particularly uncertain 
at the moment since it is even unclear, on the basis of stellar evolution theory, 
whether one expects 3 He to increase or decrease with time. On the one 
hand, 3 He has long been thought to be produced by low mass stars p8[ p9fl , 
an idea apparently supported by the observation of high 3 He abundances in 
several planetary nebulae . On the other hand, chemical evolution models 
including such 3 He production lead to overproduction of 3 He relative to the 



observations |31|, |32|, |23], [33J . 

There has been some recent progress on this question as the now venerable 
models for 3 He production have been updated [^4], jSSfl . The different groups 
agree in finding that there is net production in the standard model. It has 
been suggested [J7], however, that in fact stellar mixing might 

destroy 3 He. The inclusion of these effects certainly goes a long way in 
alleviating the problems caused by a high primordial D/H abundance p3 



f40| , fH| . Indeed, a thorough study of 3 He evolution along with other isotopes 
produced by low-mass stars would help to clarify the situation. 

Given the compounded uncertainty regarding the evolutionary history 
of 3 He from both stellar and chemical evolution, we know of no model- 
independent way to obtain a statistically significant primordial abundance 
for 3 He and therefore we do not include it in our analysis. We emphasize, 
however, that 3 He remains an important tool for chemical evolutionary mod- 
els and perhaps also for stellar evolutionary models as the destruction of 3 He 
has important consequences on the abundances of other isotopes. 



2.3 Helium-4 

The most useful site for obtaining 4 He abundances has proven to be H II re- 
gions in irregular galaxies. Such regions have low (and varying) metallicities, 
and thus are presumably more primitive than such regions in our own Galaxy. 
Because the 4 He abundance is known for regions of different metallicity one 
can trace its evolution as a function of metal content, and by extrapolating 
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to zero metallicity we can estimate the primordial abundance. 

There is a considerable amount of data on 4 He, O/H, and N/H in low 
metallicity extragalactic HII regions E3, H3I, H. In fact, there are over 



50 such regions observed with metallicities ranging from about 2-30% of 
solar metallicity. The data for 4 He vs. either O/H or N/H is certainly well 
correlated and an extensive analysis || has shown this correlation to be 
consistent with a linear relation. While individual determinations of the 4 He 
mass fraction Y have a fairly large uncertainty (AY > 0.010), the large 
number of observations lead to a statistical uncertainty that is in fact quite 
small. A recent calculation || |45| gives 



Y p = 0.234 ± 0.003(stat) ± 0.005(syst) (2) 

with similar central values and statistical errors obtained by other groups. 

We have claimed that the primordial abundance of 4 He used below is 
devoid of any serious dependence on galactic chemical evolution. The role 
chemical evolution plays here is simply to verify that the linear extrapolation 
to zero metallicity is a meaningful estimate of the primordial abundance Y p . 
The basic assumption that metal abundances increase linearly with time has 
already proven to be a secure ansatz; in addition to the strong correlation 
of the data between different metallicities, we are further assured of the 
accuracy of the estimate because the source of the data is derived from very 
low metallicity environments. 

Despite the very small statistical uncertainty in the fit shown in Eq. (|2|), 
more serious complications arise in the process of extracting abundances from 
linestrengths. As pointed out elsewhere , assumptions about the H II region 
and model-dependencies could introduce a significant systematic error. While 
the value of <7 sys = 0.005 we quote above attempts to estimate this error (see 
e.g. |4(| U), its magnitude and distribution are not well understood. Thus in 
our analysis, we will examine the effect of different assumptions concerning 
the systematic errors. 



2.4 Lithium-7 

The primordial 7 Li abundance is best determined by studies of the Li content 
in various stars as a function of metallicity (in practice, the Fe abundance). 
At near solar metallicity, the Li abundance in stars decreases with decreasing 
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metallicity, dropping to a level an order of magnitude lower in extremely 
metal poor Population II halo stars with [Fe/H] < —1.3 ([Fe/H is defined 
to be the log 10 of the ratio of Fe/H relative to the solar value for Fe/H). 
At lower values of [Fe/H], the Pop II abundance remains constant down to 
the lowest metallicities measured, (some with [Fe/H] < —3 !) and form the 
so-called "Spite plateau." With Li measured for nearly 100 such stars, the 



plateau value is well established. We use the recent results of ||7] to obtain 
the Li abundance in the plateau 

= (1.6 ±0.1) x 10- 10 (3) 

H 

where the errors are statistical only. Again, if we employ the basic chemical 
evolution conclusion that metals increase linearly with time, we may infer 
this value to be indicative of the primordial Li abundance. 

One should be aware that there are considerable systematic uncertainties 
in the plateau abundance. First there are uncertainties that arise even if 
one assumes that the present Li abundance in these stars is a faithful indi- 
cation of their initial abundance. The actual 7 Li abundance is dependent on 
the method of deriving stellar parameters such as temperature and surface 
gravity, and so a systematic error arises due to uncertainties in stellar atmo- 
sphere models needed to determine abundances. While some observers try 
to estimate these uncertainties, this is not uniformly the practice. To include 
the effect of these systematics, we will introduce the asymmetric error range 
Ai = t[j;3 which covers the range of central values for 7 Li/H, when different 
methods of data reduction are used (see eg. |f45lQ . 

Another source of systematic error in the 7 Li abundance arises due to 
uncertainty as to whether the Pop II stars actually have preserved all of 
their Li. It has been suggested that, e.g., rotational effects could reduce the 
initial Li abundance to a much lower level, and models have been advanced 
which claim to do so while maintaining the plateau behavior with respect to 



metallicity and temperature [17, pl|. While the detection of the more fragile 



isotope 6 Li in two of these stars may argue against a strong depletion P8|| , it 
is difficult to exclude depletion of the order of a factor of two. Furthermore 
there is the possibility that the primordial Li has been supplemented, by 
the time of the Pop II star's birth, by a non-primordial component arising 
from cosmic ray interactions in the early Galaxy [[|£| ^(J, 5lfl . While such 



a contribution cannot dominate, it could be at the level of tens of percent 
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4jJ . Note that this effect acts only to correct downwards the true primordial 
Li abundance from the plateau level. To allow for these possibilities, we 
will investigate the effect of a second Li systematic error, having the range 
A 2 = tglj. In fact, in § [L2| , we will turn the problem around, and by using 
the D measurements in QSOALS we will offer BBN constraints on the level 
of Li depletion. 



3 Likelihood Analysis 

Monte Carlo and likelihood analyses have proven to be useful tools in testing 
the consistency of BBN f|, |5|, [], [7|, [IUL |52| . As the observational determina- 
tions of the light element abundances improve, there is little justification for 
neglecting the uncertainties in the BBN predictions. For 4 He, these uncer- 
tainties are dominated by the uncertainty in the neutron half-life which have 
been dramatically reduced in the last several years, and lead to deviations 
(la) < 0.001. For D and 3 He, the uncertainty in the BBN yields are less 
than 10%. Of the four light element isotopes, the largest uncertainty resides 
with 7 Li, where the 1 a deviations remain as large as 20-25% 0. 

If we restrict our attention to the standard big bang model, in the context 
of the standard electroweak model with three neutrino flavors, there is one 
single unknown parameter in the standard model of big bang nucleosynthesis, 
the baryon-to-photon ratio, rj. For a given value of rj, the uncertainties in the 
BBN calculation of the light element abundances stem from the uncertainties 
in the nuclear (and weak) interaction rates employed. When we consider 
deviations from the standard electroweak theory, which can sometimes be 
parameterized by changing the number of light neutrino degrees of freedom, 
the abundances will be sensitive to the choice of N„. As we stated earlier, it 
will be sufficient to use only two of the light element abundances to test this 
one-parameter theory (one is enough to constrain it). If and when they are 
applicable, the others will even more strongly test / constrain the theory. 

The Monte Carlo calculations in BBN make available a distribution of 
abundances at each value of 77, based on the uncertainties in the nuclear and 
weak interaction rates, which we take to be Gaussian distributed. As in [|nj . 
we will use as the starting point the Monte Carlo results from Hata et al. 
[[?]]. Thus for each of the light elements we obtain a theoretical distribution 
function, which depends on 77 and on the element abundance. For example, 
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we can begin with a likelihood distribution from the BBN calculation for 
4 He: 

L BB ^Y, rj) = _U e -(>- W,)) 2 /2 CT ? (4) 



where Y BB> ^{r]) is the central value for the 4 He mass fraction produced in the 
big bang, and <J\{r\) is the uncertainty in that value derived from the Monte 
Carlo calculations. Note that there is a mild dependence on rj in o\. 

There is also a likelihood distribution based on the observations. In this 
case we have two sources of errors as discussed above, a statistical uncertainty, 
<t 2 and a systematic uncertainty, <r sys . Unfortunately, there is no well defined 
way to treat the systematic errors. One possibility is to assume that the 
systematic error is described by a top hat distribution 0, ^3|. In this case, 
the convolution of the top hat distribution and the Gaussian (to describe the 
statistical errors in the observations) results in the difference of two error 
functions 



Lq(Y,Y ( 



1 



o 



2(c"sys+ C"sys- 



erf 



Y - Yp + Qr sys _ 

V2a 2 



where in this case, Y"o is the observationally determined value for the primor- 
dial 4 He mass fraction and we have allowed for the possibility of asymmetric 
systematic uncertainty. The distribution (^) is normalized to one with re- 
spect to integration over Y. 

In addition to the top-hat distribution for the systematic uncertainty, 
we have also derived the likelihood functions assuming that the systematic 
errors are Gaussian distributed. In this case, the convolution also leads to a 
Gaussian, with an error a 2 = o\ + a 2 js . Finally, as a third case, we have also 
simply shifted the mean value Yq by an amount ±<r sys . In this case Lq is 
also a Gaussian, with spread These functions were similarly derived for 
7 Li. 

For 4 He, we constructed a total likelihood function for each value of r/ 10 = 
10 10 ?7, convolving for each the theoretical and observational distributions 

= / dYL BBN (y, rj) L Q (Y, Y Q ) (6) 
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An analogous calculation was performed for 7 Li. At this point we could use 
the two likelihood functions, and their product, to calculate (say) 68% CL 
intervals in rj. This would tell us the values of the parameter (rj) for which the 
theory (BBN) best fits the data. It tells us nothing however about whether 
that "best" fit is a good fit. For example, if the observed 7 Li abundance were 
several standard deviations below the minimum value in the 7 Li vs rj curve, 
then clearly the data fit the theory very poorly. The likelihood function for 
7 Li would still however show a peak at a value of i] near that minimum, and 
would lead to a 68% CL interval in rj similar to that obtained for a much 
higher observed value. If the only errors we had to concern ourselves with 
were Gaussian then we could calculate a x 2 from the likelihood, and this 
would then give us the probability that the theory fits the data (or more 
correctly, that the theory fails to fit the data). Unfortunately of course, the 
errors are not all Gaussian and there is no standard technique for calculating 
the goodness of fit. In this paper we estimate the goodness of fit from the 
height of likelihood functions at their peak as described below. 

We begin by defining a renormalized probability distribution in rj by a 
simple rescaling, such that 



Likelihood functions like this, derived also for lithium, form the basis of our 
subsequent analysis. We should note that, while this procedure in principle 
throws away the information on whether the observed lithium abundance is 
below the minimum value in the lithium vs rj curve (as in our example above) 
this is never an issue in practice. 

We will see that each of the abundances can individually be reconciled 
with the one-parameter theory, each predicting a distinct concordance re- 
gion in 77. Because of the independence of the observational distributions 
folded with the (dependent) theoretical distributions to give our likelihood 
distributions, these "measurements" of rj are all independent of each other. 
Demanding that two or more, in our case in particular 4 He and 7 Li, be fit 
simultaneously, constitutes a test for the theory. To do this in practice, one 
examines the product of the individual likelihoods, C = L^l^^L^^rj) 
which yields information on the goodness of fit based on the magnitude of 
this quantity at the peak of the combined distribution (if any) and of the 
spread in the allowed values in 7710. 




(7) 
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Before we move on to discuss specific results, it will be useful to discuss 
the issue of goodness of fit, and how it is to be quantified. In the case of BBN 
the number of "measurements" never increases; we will always have only 2 
to 4 of them. So there is no hope of ever being able to apply the central limit 
theorem to obtain a well defined mean and spread in the predicted range of 77 
in the usual sense of finding "world averages" of experimental quantities from 
ensembles of independent experimental results. On the contrary, in BBN one 
can always go back and refine the errors of the individual "measurements" , 
thus getting an improved test of the theory that way. Nevertheless, the 
statistical likelihood analysis will always be hampered by the small number 
of data points. Indeed, even if all of the errors are Gaussian, the usual x 2 
analysis only marginally applies for N = 2 — 4. However, in our case the 
errors and distributions are in fact far from Gaussian. Though we can still 
define a quantity \ 2 — — 21n£, its utility is unclear. However, the combined 
likelihood function does carries certain information about the goodness of fit; 
we will choose to access this by examining the significance of the amplitude. 

To introduce our approach, we first consider some simple cases to deter- 
mine what the magnitude at the peak of the combined distribution will tell 
us. For example, suppose that we we have two normalized Gaussian dis- 
tributions, Li(z) = (l/ v / 2w 1 )e- :E2 / 2<T i, and L 2 (x) = (1/ 'v^o^e - **-^ '' /2a % . 
(Think of these as toy models for L^^irf) and L j£ al (77) respectively.) The 
product of these distributions is also a Gaussian 

2 (T 2 

Li(x)L 2 (a;) = (\j2Txo x a 2 )e^ ^ 1 e^i {x ~^ ? ( 8 ) 

where a 2 = a\ + erf. In the event that the two distributions have the same 
mean value, \x = 0, then the peak of our combined distribution has a magni- 
tude l/27rcic 2 . When the distributions are offset, even though the product 
is still a Gaussian, the magnitude at the peak is now suppressed by a factor 
e -/i /2cr _ Clearly for [i 3> a, this suppression can be significant and indicates 
a lack of goodness of fit. It would seem appropriate, therefore, to compare 
the value of the peak of the combined likelihood distribution (^j) with the 
quantity (l/27r<7i<72). Despite the fact that the suppression factor in @ is 
exponential, unless // ^> a, it will take values in the tens of percent. For 
example, if o\ ~ cr 2 and if our distributions are offset by 1 (2) a, the suppres- 
sion factor is still relatively mild, 0.78 (0.37). Even for a distribution whose 
peaks are separated by 3a, i.e. distributions that we would generally judge 
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as being inconsistent, the suppression is only about 0.1. Thus when we make 
the comparison at the peak of our combined distribution to (l/27r<7i<72), any 
value significantly less than 1 will indicate a poor fit. 

Thus far in our example, we have only considered the case in which the 
likelihood functions, L(r]) are Gaussian. However, in general our likelihood 
distributions for 4 He and 7 Li will not be Gaussians. In particular, as we will 
see, in the case of 7 Li, our distribution is double peaked. We can however 
still make the comparison to (l/27ro"iO"2) if we define the respective cr's as the 
half-width at half-maximum of a local peak. This will be made clear when 
we consider specific examples below. 

In considering the bimodality of the Li likelihood function, it is useful 
to adopt a larger perspective. Over an extended range of 77, beyond the 
canonical 7710 ~ 1 — 10 range,[] all of the light element abundances are non- 
monotonic functions of r]. Thus, including a larger range of r\ in one's analysis 
would yield multiple peaks in each isotope's individual likelihood function. 
Of course, the "new" peaks would not overlap between different isotopes, and 
the predicted range for rj as given by the combined likelihood function would 
remain unchanged. Hence it is reasonable to independently normalize each 
well-isolated peak in an isotopic likelihood. This procedure is straightforward 
for all of the isotopes other than Li, whose peaks are not isolated and can 
overlap to some extent. Thus we normalize over both of these peaks and 
therefore expect a reduced amplitude in the combined distribution. 

4 Model-Independent Odds on BBN 

4.1 Standard Cases 

The procedure that we set up in the preceding section would be straightfor- 
ward to carry out if it were not for the complication of the treatment of the 
systematic uncertainties in the data. We will therefore present results for 
various different assumptions regarding the data. As described in section 2, 
we will always make comparisons with respect to our standard case in which 
the data is described by Eqs. (|) and (|), namely Y p = 0.234 ± 0.003 ± 0.005 

1 Indeed, if one is testing BBN itself, one should in principle consider a large range in 
77, but of course experience with BBN calculations, as well as other available bounds on 
f?B and hence on 77, lead one to focus straightaway on the canonical range. 
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and 7 Li/H = (1.6 ± 0.1 ±i) x 10~ 10 . The effects of lithium depletion are 
neglected at this stage, and the systematic errors are described by a top-hat 
distribution. 

With these assumptions, we have calculated the likelihood functions for 
4 He and 7 Li (TTJ] which are shown in figure [|. The shapes of these curves are 



characteristic, with one peak for 4 He (whose abundance rises monotonically 
with rj), and two for 7 Li (whose abundance goes through a minimum). In 
this case (and most others) the minimum theoretical Li is somewhat below 
most of the observational values and so the sides of the minimum are favored, 
leading to the two peaks; i.e., for a given observational value of 7 Li, there are 
two values for rj at which this may be achieved. By glancing at figure [l], it does 
not take very much statistical machinery to see that the BBN predictions for 
4 He is consistent with the low-77 value of the 7 Li prediction when compared 
with the observations. These distributions are clearly consistent. 

The combined likelihood, for fitting both elements simultaneously, is given 
by the product of the two functions in figure [l], and is shown in figure j|. As 
we discussed in the previous section, we have scaled the combined distribu- 
tion by a factor 27r<74<77 where a 4 = .66 is the half-width at half-maximum 
of -^totail 7 /) an d 07 = .38 is the corresponding quantity for L^^rj). In this 
latter case, the distribution i^totai^) * s clearly far from Gaussian and the 
value chosen for cr 7 corresponds to the low-77 peak only. Because the 7 Li 
distribution is nearly equally divided into two peaks we should expect com- 
pletely overlapping distributions to yield a peak value of only ~ 0.5 in the 
combined distribution^]. This is what one finds in figure which does show 
concordance, the peak is indeed close to 0.5 at r] 10 = 1.8. 

The allowed 68% CL and 95% CL ranges are 

1.6 < 7710 < 2.8 68%CL 

1.4 < J7 10 < 3.8 95%CL (9) 
Note that these intervals are based on standard statistical techniques, and 



2 Note that the resulting peak value of 0.5 does carry some statistical information. For 
example, suppose our 7 Li distribution takes a form in which the peaks were different from 
one another (e.g. if the left peak were a factor of 10 smaller than the right peak). While 
we could argue that we would expect a combined fit of only 0.1, we should also conclude 
a low confidence level for consistency. Thus it is not appropriate to scale out this factor 
in the combined distribution function. However, the Li peaks are nearly always about 
equal. 
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are not dependent on our chosen method of assigning a goodness-of-fit to the 
height of the likelihood peak. 

Thus, for this "standard" case, we find that the abundances of 4 He and 
7 Li are consistent, and select an rji range which overlaps with (at the 95% 
CL) the longstanding favorite range around 7710 = 3. Further, by finding 
concordance (in this case) using only 4 He and 7 Li, we infer that if there 
is problem with BBN analysis, it must arise from D and 3 He and is thus 
tied to chemical evolution. The most model-independent conclusion is that 
standard BBN with N v = 3 is not in jeopardy, but there are problems with 
our detailed understanding of D and particularly 3 He chemical evolution. 

The concordance range for 77 given in Eq.(|D allows us to make definite 
predictions for the primordial abundances of D and 3 He. The 95 (68) % CL 
ranges in © corresponds to 5.5(8.9) < (D/H) x 10 5 < 28(22) with a best 
value for D/H = 1.8 x 10~ 4 at i] w = 1.8. As we have already noted, this 
value for D/H agrees incredibly well with the QSOALS observations in [50 . 
For 3 He, we have, 1.4(1.7) < ( 3 He/H) x 10 5 < 2.7(2.5), with a best value 
3 He/H = 2.3 x 10~ 5 which, it should be noted, is larger than the solar value 
of 3 He/H. 

In the preceding analysis, we have described the systematic uncertainty 
by a top-hat distribution. This prescription gives a tight set of errors which 
provides a stringent test to BBN, however, it also gives the most restrictive 
bounds on rj. Thus when we relax this assumption and treat the systematics 
as Gaussian distributed, we will expect concordance over a greater range. In 
figure |3], we show the 4 He and 7 Li likelihood functions when the systematic 
errors are Gaussian. It should not be surprising that the concordance is 
present at the same level though the widths of the distributions are larger. In 
this case, 04 = .86 and 07 = .53. The combined distribution shown in figure § 
is also similar to that in figure |2|, the peak is again close to 0.5 and the width of 
the distribution is larger, placing greater weight at higher 77. This distribution 
leads to a predicted range for r], 1.4(1.7) < r] w < 4.3(3.4) The peak of the 
distribution is unchanged at rj w = 1.8. This range in 77 corresponds to 
4.6(6.5) < D/H x 10 5 < 28(20) and 1.3(1.5) < ( 3 He/H) x 10 5 < 2.7(2.4). 

Having found concordance for BBN in the standard model, we now vary 
the parameters of the distributions — in particular, those for the systematic 
errors — and examine the results. 
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4.2 7 Li Depletion 

Next we would like to test the uncertainty in 7 Li given by A 2 . That is, what 
is the effect of 7 Li depletion and cosmic-ray production? We will first assume 
that the 7 Li/H abundance is shifted up by A 2 , that is, we assume a factor 
of 2 depletion in the Pop II halo stars. We will treat the error in Ai, as well 
as the systematic error in 4 He as Gaussian, this being a more conservative 
approach than the top-hat method. 

By shifting the observational primordial 7 Li abundance upwards, we are 
in effect separating the two peaks in the 7 Li likelihood distribution. The 
results for this case is shown in figure [f| There is still some overlap between 
the 4 He and 7 Li distributions. In this case however, there is a distinction 
between the two lithium peaks. The half- width of the low-?? peak is 0.19 
while for the high-?] peak it is 0.68. In the combined distribution, we have 
scaled the overlap regions according to these respective peaks, which can be 
justified in this case since the combined distribution now shows two distinct 
distributions as seen in figure []. This combined likelihood distribution is in 
clear contrast to the previous cases considered. Note the peak value at low-77 
is only about 0.15 (it is even lower for the high-?/ peak). We do not consider 
this to be a good fit of theory to data; yet on the basis of this result relying 
4 He and 7 Li alone, we cannot quite exclude depletion at this level. 

To salvage the consistency of BBN with a factor of 2 depletion in Li one 
might consider an increase the systematic uncertainty in 4 He by a factor of 
2, to cr sys = 0.01. By doing so, we would flatten and greatly broaden the 
4 He likelihood distribution. This would of course produce a greater overlap 
with the high-?/ peak in the 7 Li distribution. However, we can not place any 
great significance to the improved overlap in this case, because by increasing 
the systematic uncertainty to this extent, we have in effect taken all of the 
predictive power away from 4 He. ( The 4 He distribution in this case is now 
mildly peaked at rj 10 ~ 1.8 at a value of 0.2, and is only slowly decreasing 
out past 7710 = 10. ) Thus we are only really considering 7 Li in this case, and 
some level of consistency is almost trivially guaranteed. 

Another possibility to salvage BBN consistency with 7 Li depletion would 
be to shift up the primordial 4 He abundance. However, a shift up by an 
amount <r sys = 0.005 to a primordial value Y p = 0.239 ± 0.003 results in a 
terrible fit in which the combined distribution peaks at a value of 0.05! Only 
if the primordial abundance is shifted by at least 0.01 is the overlap between 
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4 He and 7 Li significant enough to give a reasonable combined distribution. 
It is our opinion, however, that in this case we are straying far from the raw 
observations which show the remarkably good compatibility in figures 1-4. 

Let us also comment that if one assumes that perhaps 20% of the observed 
7 Li were produced by cosmic-ray nucleosynthesis rather than BBN, then 
compatibility is further improved. If we shift the 7 Li down by A 2 , the two 
7 Li peaks begin to merge and there is a strong overlap in the distribution 
functions. The combined distribution peaks at 7710 = 2.0 in this case. 

Finally, we point out that in order to bring the probability of overlap of 
the theoretical and observational values of the Li abundance]] below the 5 per 
cent level, would require reducing the observationally inferred abundance by 
an order of magnitude. This would be the lowest observational value BBN 
could be in agreement with, but at present the limit does not seem at all 
interesting. 



4.3 Primordial D/H 

It is interesting to note that the central (and strongly) peaked value of rji 
determined from the combined 4 He and 7 Li likelihoods is at r/ 10 = 1.8. The 
corresponding value of D/H is 1.8 xlO -4 , very close to the value of D/H in 



quasar absorbers in the published set of observations jlTJ, It is not clear 
if this is a coincidence or if we really have evidence that three of the light 
element abundances point to the same value of 7710. 

Noting the perhaps still preliminary nature of the QSOALS D/H mea- 
surements, we move ahead and perform the likelihood analysis for the three 
light elements D, 4 He and 7 Li. To include D/H, we proceed in much the same 
way as with the other two light elements. We compute likelihood functions 
for the BBN predictions as in Eq. (^) and the likelihood function for the 
observations using D/H = (1.9 ± 0.4) x 10~ 4 . In this case, since the system- 
atic error is not stated, it is neglected here and we adopt a simple Gaussian 
form for the statistical errors. These are then convolved as in Eq. (|6|). In 
figure 7, the resulting normalized distribution, Lf otal (rj) is super- imposed on 
distributions appearing in figure 1. It is indeed startling how the three peaks, 
for D, 4 He and 7 Li are literally on top of each other. In figure 8, we show the 

3 This is given by an integral equation similar to (7) for lithium with the lower limit of 
integration set to the theoretical minimum and letting the observational value vary until 
the value of the integral becomes 0.05. 



19 



combined distribution which has now been scaled by ((27r) 3//2 <720"4<77) 1 with 
o"2 = -29. We have a very clean distribution and prediction for rj: 

1.65 < 7710 < 2.06 68%CL 

1.50 < ?? 10 < 2.37 95%CL (10) 

with the peak of the distribution at 7710 = 1.75. The absence of any overlap 
with the high-7? peak of the 7 Li distribution has considerably lowered the 
upper limit to 77. Overall, the concordance limits in this case are dominated 
by the deuterium likelihood function, which, we should caution again, is 
based on an observation along a single line of sight. 

A high value of D/H, as measured in the QSOALS, requires a signifi- 
cant amount of D destruction over the history of the Galaxy. Although this 
alone is not necessarily problematic, since chemical evolution models can 



be constructed to account for such factors of deuterium destruction |31| , it 
compounds the problem of 3 He overproduction. The fact that 3 He was prob- 
lematic even at higher r] w ~ 3, was our motivation for neglecting 3 He to 
begin with. 

Finally, we consider the effects of the D/H measurements when 7 Li is 
assumed to be depleted. In figures 9 and 10, we show the individual and 
combined likelihood functions which can be compared with those in figures 5 
and 6 with D/H. In the combined distribution, the high-?] peak is now gone, 
and the low-?? peak is hardly significant. Recall (§ [4.2| ) that 7 Li depletion could 
be tolerated by either an enlargement of the 4 He systematic error or a shift 
upwards in the 4 He abundance, which would would allow significant overlap 
with the high-77 peak. However, when the D/H observations are included in 
the analysis, the high-?? peak is absent in the combined distribution and the 
low-?? peak takes a value of less than 0.05. Once the D/H measurements are 
confirmed, they will provide a strong constraint on the degree of 7 Li depletion 
in halo stars when the consistency of big bang nucleosynthesis is assumed. 



5 Constraints on N u 

In the previous section we have demonstrated the health of BBN when model- 
independent constraints are applied; we arrived at these conclusions without 
first assuming the validity of BBN. Now we assume its validity and, as has 
been done traditionally, use it to constrain new physics. In particular, we 
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will consider the effects of our analyses on the limit to the number of neutrino 
flavors, iV,,. 

The light element isotope which is most sensitive to the number of neu- 
trino flavors is 4 He. In fact, it has been common [l[] to express (by means of a 
fit) the 4 He mass fraction as a function of r] w , N v , and the neutron meanlife. 
For example, near r] w = 2 and N v = 3, we have, 

Y = 0.2262 + 0.0131(^ - 3) + 0.0135 lnr/ 10 + 2 x 10" 4 (x„ - 887) (11) 

Therefore, given an upper limit to the 4 He mass fraction and a lower limit 
to r\ one can derive an upper bound to N„. Indeed, unless the lower bound 
to t]io is greater than 0.4 (something we are assured of now by the D/H 



measurements), there is no upper limit to N u [p4|] . 111 H3; h was argued 
that the solar value of D+ 3 He could provide a reasonable lower bound to 
r\. Recall that this argument was made before the evidence from planetary 
nebulae pointed towards the need for 3 He production in low mass stars and 
before the QSOALS measurements of extremely high D/H. Therefore, the 
role of chemical evolution was not believed to be critical. Upper limits to N v 
were then derived by inserting the 2 a upper limit to Y, the lower bound on 
7] from D + 3 He, and the lower bound on r n into Eq. (|Tl~|). The result found 
in was N u < 3.3. 

The methods for deriving N„ have since been refined incorporating a more 
statistical meaning to the bound. In [||, using the best central values and 
their associated uncertainties in an expression similar to ([11]) gave a best 
estimate, N u = 2.2 ± 0.3 ± 0.4, (where the statistical error is taken from 
the observational determination of Y and the neutron mean life, and the 
systematic error from 4 He and from t], assuming that rjio = 3.0 ± 0.3). Since 
one could well imagine theories that would effectively lower the value of N v as 
well as increase it, one has to, in the broadest sense take the bound seriously, 
and accept that it might show preference for some extension of the standard 
model predicting less helium. The (2 o" sta t + o~ sys ) upper limit was found to 
be N v < 3.1 PJ. A similar bound, N v < 3.04 was obtained from a Monte 
Carlo analysis |J and an even stronger bound was obtained using a more 
refined treatment of the solar D + 3 He argument, N v = 2.0 ± 0.3 [^2]| . If we 
assumed a priori that N v > 3, say, then we could use a Bayesian approach to 
relax the upper limit to N v [[53| , p5[| . An alternative argument for relaxing the 
constraint based on a stochastic approach to chemical evolution was made 
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in | 5"E|| . Also, for a critical review on the history of the bound on N v see also 



ES- 
AU of the above methods for placing a bound on N u rely on 3 He and in 

particular, the D + 3 He argument for a lower bound to rj. However, as we have 

emphasized repeatedly in this work, the uncertainties in the evolutionary 

history of 3 He make it a poor indicator of N u or of the validity of BBN. If 

instead we use the results found in §4, for rj based on 4 He and 7 Li, we find 

as the best estimate for N v , 

N v = 2.99 ± 0.23 ± 0.38 (12) 

showing no particular preference to N u < 3, in fact preferring the standard 
model result of N u = 3 and leading to N u < 3.90 at the 95 % CL level 
when adding the errors in quadrature. (Since the central value in Eq. (|12]) 
is dramatically close to 3, we show two decimal places, but of course the 
errors imply that the precision of the agreement is somewhat fortuitous.) In 
(|T2"D, the first set of errors are the statistical uncertainties primarily from the 
observational determination of Y and the measured error in the neutrino half 
life T n . The second set of errors is the systematic uncertainty arising solely 
from 4 He, and the last set of errors from the uncertainty in the value of rj and 
is determined by the combined likelihood functions of 4 He and 7 Li, ie taken 
from Eq. (|9|). We view Eq. ( |T2| ) as a further (and remarkable) confirmation 
of standard BBN. 

Since we have shown our likelihood results when quasar D/H is included, 
we can also derive the best value for N v when D, 4 He and 7 Li are all included 
in the analyses. In this case, 

N v = 3.02 ± 0.23 ± 0.38 1^ (13) 

leading again to a 95 % CL limit of N v < 3.9. In Eq. (|I3"1) , the central value is 
slightly higher due to the central value for 7] w being slightly lower in this case 
(1.75 rather than 1.80) and the error due to the uncertainty is considerably 
smaller especially on the low side (corresponding to the high side on 77). 



6 Conclusions 

We have argued that because D and 3 He are observed only in the solar sys- 
tem or the local ISM, determinations of their primordial abundances are 
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particularly uncertain and dependent on models of galactic chemical evo- 
lution. It is therefore important to determine the status of BBN without 
them. Using a likelihood analysis we find that the compatibility of the 4 He 
and 7 Li constraints depends on the assumptions one makes about the size 
and distribution of the possible systematic errors. Nevertheless, for the most 
"standard" assumptions, i.e. taking the data on 4 He and 7 Li at face value, we 
find that these two isotopes are quite compatible and give BBN concordance 
to a high degree of confidence. The baryon-to-photon ratio selected by this 
concordance is centered around r] W ~ 1.8 and ranges up to 3.8 (or 4.3 if all 
errors are assumed to be Gaussian). 

While the putative determinations of D in quasar absorption line systems 
are still in their infancy, the technique can be a powerful probe. Upon intro- 
ducing into our analysis the relatively high D abundance reported for several 
systems, we find that the agreement with 4 He and 7 Li is striking. It is crucial 
to further observe and better understand these systems, which can provide 
a decisive constraint on BBN. 

Having found BBN theory to be consistent with observation, we turn to its 
implications. In the realm of astrophysics, we find that the needed primordial 
D and 3 He abundances challenge our understanding of the chemical evolution 
of these nuclides, particularly 3 He. We find that the presence of 7 Li depletion 
does not improve the BBN agreement, and further, if one adopts the high 
D abundance suggested by the quasar measurements, one can determine the 
primordial 7 Li abundance and thus tightly constrain the degree of depletion. 
In the realm of particle physics, we revisit the N u counting argument, with a 
careful treatment of all the data giving a best value for N v = 3.0 and a 95% 
CL upper limit of N v < 3.9. 
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FIGURE CAPTIONS 



1. Normalized likelihood distribution for each of 4 He and 7 Li, shown as a 
function of rj. The one-peak structure of the 4 He curve corresponds to 
its monotonic increase with 77, while the two-peaks for 7 Li arise from its 
passing through a minimum in the theoretical calculation. Systematic 
uncertainties have been assumed to be top-hat distributions. 

2. Combined likelihood for simultaneously fitting 4 He and 7 Li, as a func- 
tion of r). Plotted is the product of the normalized likelihood distribu- 
tions shown in figure [3] multiplied by 27r0" 4 cr 7 . 

3. As is Figure 1, where the systematics have been assumed to be Gaussian 
distributions. 

4. As in figure 2, using the likelihood distributions in figure 3. 

5. As is Figure 3, where a factor of 2 depletion for 7 Li has been assumed. 

6. As in figure 2, using the likelihood distributions in figure 5. 

7. As in figure 1, adding the likelihood for D/H as given by the QSOALS 



observations IS 



8. As in figure 2, using the likelihood distribution in figure 7. 

9. As in figure 5, assuming the observed Pop II 7 Li to be depleted by a 
factor of 2, and including the D/H as observed in QSOALS. 

10. As in figure 6, using the likelihood distribution of figure 9. 
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